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Abstract 



Context. The physical mechanisms driving angular momentum transport in accretion discs are still unknown. Although it is generally 
accepted that, in hot discs, the turbulence triggered by the magneto-rotational instability is at the origin of the accretion process in 
Keplerian discs, it has been found that the values of the stress-to-pressure ratio (the a "viscosity" parameter) deduced from observa- 
tions of outbursting discs are an order of magnitude higher than those obtained in numerical simulations. 

Aims. We test the conclusion about the observation-deduced value of a using a new set of data and comparing the results with model 
outbursts. 

Methods. We analyse a set of observations of dwarf-nova and AM CVn star outbursts and from the measured decay times determine 
the hot-disc viscosity parameter ah. We determine if and how this method is model dependent. From the dwarf-nova disc instability 
model we determine an amplitude vs recurrence-time relation and compare it to the empirical Kukarkin-Parenago relation between 
the same, but observed, quantities. 

Results. We found that all methods we tried, including the one based on the amplitude vs recurrence-time relation, imply a h ~ 0. 1 -0.2 
and exclude values an order of magnitude lower. 

Conclusions. The serious discrepancy between the observed and the MRI-calculated values of the accretion disc viscosity parameter 
a is therefore real since there can be no doubt about the validity of the values deduced from observations of disc outbursts. 

Key words. Accretion, accretion disks - Instabilities - Stars: dwarf novae 



1. Introduction 

Since the very beginning of the accretion-disc theory, the mech- 
anism of angular momentum transport through the disc has been 
a matter of debate, and despite strenuous efforts of many re- 
searchers it remains an open problem because the results of nu- 
merical simulations do not match observati ons when the respec- 
tive v alues of the vi scosity paramete r a (Sha kura & Sunvaevl 
1 1973b are compared dKing et al.l l2007, ; see, however, Sorathia 
et al. 2012). Magneto -rotational instability (hereafter MRI; 
Balbus & Hawlevl fl998) simulations result in average a values 
of the order of 0.01, whereas the best studied case of dwarf- 
nova outbursts unambiguo usly pr ovide values that are an order 
of magnitude higher dSmakl fl999, hereafter S99). 

In addition, from simulations of dwarf nova (DN) eruptions 
it is clear that a must vary during outbursts. The outburst am- 
plitudes can be reproduced only if a c in the quiescent disc is 
four to ten times sma ller than g;, in the hot, outbursting state 
(see e.g.[Smak 1984at iMever & Mever-Hofmeisted[T984 . This 
holds regardless of the disc chemical composition (however, the 
required a jump i s rather 2-6 for helium-dominated discs, see 
iKotko et alj2012l) . Since observations impose as 0.1-0.3, the 
cold disc a c should be ~ 0.01. This is incons istent with the MRI 
simulations that result in gh ~ 0-01 (see e.g.lHirose etal.l [2009) 
and do not predict a jumps dSorathia et alj|2012l obtain a jumps 
but with timescales too short to correspond to the dwarf-nova 
case.) In the present article we revisit the problem of deter- 
mination using a new set of data and methods that are somewhat 
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different from those used by S99. We confirm his general con- 
clusions about the value of g h - 

The values of the viscosity parameter a deduced from obser- 
vations of dwarf-nova decay from outburst maximum are almost 
model-independent. Basically one only assumes that the decay 
time reflects the viscous character of this process. 

On the other hand, the phenomenological relation be- 
tween the outburst amplitude and the recurrence time 
(the so-called Kukarkin-P a renago relatio n, hereafter K-P; 
iKukarkin & Parenagol fl934t I Warned f2003) obviously reflects 
some properties of the outbursts and should be derivable from 
the model supposed to be describing dwarf-nova outburst s, i.e. 
from the disc instability model (hereafter DIM; see lLasotal 
1200 ll for a review). Using the DIM we derive an amplitude - 
recurrence-time relation that compares reasonably well with the 
K-P relation, especially considering the large scatter of the ob- 
servational data. We find that this K-P type relation also implies 
that g h cannot be of the order of 0.01 but must be roughly ten 
times larger than this value. 

In Section [2] we briefly describe those aspects of the DIM 
that are required for understanding dwarf nova outbursts and the 
quantities that characterize them. The estimate of g h from the 
empirical relation and from the analytical derivation of the out- 
burst decay time are presented in Section[3] Section|4]deals with 
the details of the analytical derivation of the Kukarkin-Parenago 
relation and its consequences for the parameter g h . Conclusions 
are presented in Section[5] 
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2. The dwarf nova outburst cycle in the DIM 
description 

Outbursts of DN can be divided into two types: normal out- 
bursts and superoutbursts. As "normal" one defines outbursts 
that are narrow, i.e. lasting no longer than about ten days, and 
that have amplitudes lower than ~ 6 mag and during which no 
superhumps (low-amplitude brightness variations) have been ob- 
served. Superhumps are the identification mark of the superout- 
bursts, which are also brighter and last longer than normal out- 
bursts. In what follows "outbursts" only refers to normal DN 
outbursts. DN showing only normal outbursts are classified as 
U Gem-type, whereas binaries exhibiting both normal and super 
outbursts belong to the SU UMa type. 

In many DN of U Gem-type (e.g. in SS Cyg and U Gem), a 
bimodality in the normal outbursts widths has been observed. 
The "wide" outbursts are ~ 0.2 mag brighter and evidently 
longer than the "narrow" outbursts. It has not been defined how 
much longer the outburst should be to be recognized as "wide", 
but the difference is clear when looking at a specific light-curve. 
The wide outbursts should not be confused with superoutbursts, 
the former being normal DN outbursts with no superhumps de- 
tected, which in turn are the characteristic feature of the latter. 
According to the DIM, wide outbursts correspond to outside-in 
heating front propagation, while narrow outburs ts resu lt from the 
opposite sense of motion of these fronts dSmak« 1984a). 

During quiescence a dwarf nova disc accumulates matter un- 
til somewhere the temperature crosses the value critical to the 
onset of the thermal instability. A heating front starts propa- 
gating into the low-temperature regions, leaving behind ionized 
matter. If the radius where it started is close to the inner disc edge 
the front will propagate in the outward direction and the outburst 
will be of the inside-out type. If the front travels inwards from 
the outer edge, the outburst is an outside-in outburst. 

In the post-front hot regions of the disc the angular mo- 
mentum transport is defined by and the accretion rate M accr 
rises. The enhanced efficiency of the outward angular momen- 
tum transport in the hot disc (ah>a c ) causes mass, which during 
the quiescence had gathered in the outer parts of the disc, to dif- 
fuse inwards at a high rate changing the surface-density 2 profile 
from 2 ~ R to 2 ~ R~ 3 I 4 . Once the heating front arrives at the 
outer (inner) disc edge, matter across the whole disc becomes 
ionized, and M accr everywhere in the disc becomes roughly con- 
stant and higher than the mass-transfer rate from the secondary 
(See Fig. @J. In thi s phase the outburst is at its maximum, the 
disc is hot and the lShakura & Sunvaevl (1 19731) solution is a good 
approximation of such a quasi-stationary configuration. 

Because of ongoing accretion onto the central body the ac- 
cretion rate (and density) decrease at a characteristic viscous 
timescale. Because the hot-disc critical surface density 2^ (/?) ~ 
R • (see AppendixlAl). 2 on the outer edge easily falls below its 
critical value and a cooling front starts propagating inwards (see 
Fig. [3J. The decay from outburst maximum can be seen as due 
to viscous emptying of a hot disc with a shrinking outer radius. 

Once the disc has lost matter accumulated during quiescence 
the refilling process restarts. Strictly speaking the disc is filling 
up already during the outburst what allows the existence of out- 
burst cycles with practically no quiescence. 

3. The value of the viscosity a parameter in the hot 
accretion discs 

As mentioned above the decay timescale is determined by the 
viscous timescale f v ; sc in the hot disc. Since f v ; sc ~ R 2 /v, assum- 



ing that the kinematical viscosity coefficient v ~ ac s H, where c s 
is the sound speed and H the disc semi-thickness, determining 
the decay time allows estimation of a, assuming one can esti- 
mate the disc radius (e.g. from the binary's orbital period). 

S99 obtained an observational relation between the outburst 
decay rate Td ec (or alternatively the outburst width W) and the or- 
bital period P or b: fdec^orb) and W(P or b). After comparing it with 
the relations of the same type found from the fits to data from 
the numerical models, he concluded that the best agreement be- 
tween observation and models is obtained for a^- 0.2. Below, 
we determine the viscosity parameter using a different DN sam- 
ple and different numerical models. 

3.1. a from the decay rate vs. the disc radius relation 

To be certain that the result is independent of the choice of the 
DN sample for deriving the Td ec - Pmb relation, we used a set of 
observationa l data d iffer ent from S99, who used the data from 
Ivan Parad iis (1983) and IWarnerl (120031) that include all types 
of DN: the U Gem, S U UMa, Z Cam stars, while our sample 
contains systems from lAk et al.l d2002l) . together with an addi- 
tional seven U Gem-type DN for which we measured the decay 
rates using light-curves from the AFOEV database. Among the 
twenty one systems that we took into account nine have been 
included also in S99. 

Since it has been shown that the basic DIM (with no mass 
transfer rate enhancement or the additional source s of disc heat- 
ing) is only able to reproduce normal outbursts (see lLasotafe OOl: 
iKotko et al.ll20H and references therein), we find it reasonable 
to start our analysis with U Gem-type systems where only nor- 
mal outbursts are presenQ to ensure that the comparison between 
Tdec measured in the models and in the real light-curves is con- 
sistently defined. Then, to test to what extent other DN types 
and the disc chemical composition influence the correlation be- 
tween Tdec and P or b, normal outbursts of SU UMa stars and one 
AM CVn system (PTF1J0719) have been added to the sample. 

The systems considered are listed in Table Q] and marked in 
Fig. [TJ The decay rates of U Gem-type binaries (rows 1-8) and 
the decay rates for the normal outbursts in SU UMa-type binaries 
(rows 14 - 20) are taken from lAk et al.l(l2002l). In rows 9 - 1 3 are 
7 U Gem-type binaries not included in lAk et aLl £2002). For each 
of them we have measured Td ec as the time it takes the system 
brightness to decline by ~ 1 mag starting from the level 1 mag 
below the maximum. The last row stands for the only AM CVn- 
type star for w hich the existence o f the normal outburst have 
been confirmed dLevitan et alJl201 Th . 

First, we compared our fit to the Td ec - ^orb relation with 
the one obtained by S99. The linear fit to our data of the form 
Tdec = C T P or b (the "Bailey relation"; Bailevl [l975l) gives C T = 
0.37 + 0.03 with dispersion rms = 0.698 which compares nicely 
with the result of S99: C T;S99 = 0.38 ± 0.02 with rms = 0.54. 
For the more general case when Td ec = C l T P ?f we get C\ = 
0.69 + 0.17 and p = 0.66 + 0.14 with dispersion rms = 0.64, 
while S99 result is C\ S99 = 0.61 ± 0.07 and p S99 = 0.71 ± 0.09 
with dispersion rms = 0.48. Our results are similar to those of 
S99 within the uncertainty of the fit coefficients. 

Next, we compared model decay times with those observed 
during dwarf-nova outbursts. Because models are calculated not 
for a given P or \, but for a given disc radius, we convert orbital 
periods of observed systems into the disc radii. 



1 In U Gem itself, one superoutbur st has been observed (e.g. 
Maso n et al . 1 988 ; Smak & Waaeen 2004) 
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Table 1. Decay-from-outburst properties of selected dwarf novae. 





System 


^orb 


Mi 


M 2 




W 


Tdec 


TdecQ 2 






(hr) 


(M ) 


(Mo) 


(10 10 cm) 


(d) 


(d/mag) 


(d/mag) 


1 


BV Cen 


14.67 


1.24 


1.1 


9.76 


20.9 


5.7 


5.47 


2 


AT Ara 


9.01 


0.53 


0.42 


5.354 


4.1 


2.3 


1.98 


3 


RUPeg 


8.99 


1.21 


0.94 


7.05 


7.2 


3.2 


2.75 


4 


MU Cen 


8.21 


1.2 


0.99 


6.586 


7.9 


3.1 


4.15 


5 


SS Cyg 


6.6 


0.81 


0.55 


5.069 


6.4 


2.5 


1.89 


6 


TWVir 


4.38 


0.91 


0.4 


4.168 


3.7 


0.8 


1.82 


7 


SS Aur 


4.33 


1.08 


0.39 


4.462 


4.3 


1.6 


1.48 


8 


U Gem 


4.25 


1.2 


0.42 


4.58 


3.2 


1.3 


1.32 


9 


EY Cyg 


11.02 


1.1 


0.49 


8.201 


14.3 


4.29 


2.83 


10 


DX And 


10.57 


1.2 


0.8 


7.923 


12.5 


3.47 


2.90 


11 


EXDra 


5.04 


0.75 


0.56 


4.097 


? 


2.63 


2.01 


12 


BDPav 


4.30 


1.15 


0.73 


4.307 


? 


2.96 


1.54 


13 


IP Peg 


3.797 


1.16 


0.55 


4.077 


? 


2.16 


1.75 


14 


CU Vel 


1.884 


1.23 


0.15 


2.988 


4.5 


1.1 


1.05 


15 


WXHyi 


1.796 


0.9 


0.16 


2.519 


4.0 


0.9 


1.55 


16 


ZCha 


1.788 


0.84 


0.13 


2.503 


3.7 


1.0 


1.24 


17 


VW Hyi 


1.783 


0.67 


0.11 


2.292 


4.0 


0.7 


0.90 


18 


OYCar 


1.515 


0.64 


0.086 


2.07 


4.7 


0.8 


0.74 


19 


EkTrA 


1.509 


0.46 


0.09 


1.775 


3.0 


0.7 


0.50 


20 


SW UMa 


1.364 


0.71 


0.1 


1.987 


7.6 


0.6 


0.49 


21 


PTF1J0719 


0.446 


0.5 


0.05 


0.871 


1 


0.25 


0.79 



Notes. P orb is the orbital period, M t and M 2 the primary and secondary masses respectively, R^ mali the maximum disc radius, W the observed 
outbur st width, Td ec the observed decay rate, and Tdecoo the theoretical (Eq.[7} decay rate for = 0.2. Rows: 1 - 8: U Gem stars from A k et al.l 
( 120021) . 9- 13: U Gem s tars with T dec c alculated by authors from light-curves fr om AFOEV da t abase 14-20: SU UMa stars from lAk et alj 12002), 
21: AM CVn star from lLevitan et al.l d201 lh : P OI b, M\ and M 2 are taken from lRitter & Kolbl d2003h : ifdmax is calculated according to Eqs. Q} & 



During outburst the outer disc ra dii in DN expand up to 
^Amax ~ 0.9Rl 1 (see e.g. Smak 2001), where R]., is the radius 
of a primary Roche-lobe given by the Eggl etonl ( 1 19831) formula: 

R U _ 0.49 g 1 ^ 

a 0.6 ? 2 / 3 +log(l + q 1 ^) K ' 

a = 3.5 x 10 ll) M2 /3 (l + q) l/3 P 2 ^ 3 cm (2) 

where q = Mx/Mo and Phi is the orbital period in hours. 

The primary and the secondary masses (M\ and M2) for cal- 
culations of Rd max fr om P nr h were taken f rom the latest version 
(November 201 1) of iRitter & Kolbl J2003h . For PTF1J0719, M x 
and M2 have been guessed according to what is expected for 
AM CVn stars since no observational estimates have been sug- 
gested yet. The fit in the form T dec = ^l^d.max to all data from 
Table[TJgives A\ = 0.48 + 0.02 with dispersion rms = 0.56. 

It is interesting to check whether the T de c - Pd.max relation is 
independent of the class of systems exhibiting normal outbursts. 
The T"de C -Pd,max relation for U Gem-type binaries is linear to a 
good approximation. For the linear fit one obtains A\ = 0.49 + 
0.03 with dispersion rms = 0.68. The coefficients for the general 
relation Td ec = B\R^ m!UL y are B\ = 0.27 ± 0.12 and y = 1.32 + 
0.22 with rms = 0.66. In this case, similar to what was noticed 
by S99, the rms dispersions do not differ significantly for the 
linear and non-linear fits, moreover the errors of B\ and y are 



rather large. Therefore we limit our further considerations to the 
simpler linear case. 

According to the model, the outbursts appearing between su- 
peroutbursts of SU UMa stars have the same origin as those in 
U Gem- type binaries. As expected, their measured Td ec marked 
on the Tdec-Pd,max plane extrapolate the Td ec -Pd ; max relation for 
U Gem to the regime of orbital periods shorter than 2 hr. The 
coefficient A\ of the linear fit for the sample, including the nor- 
mal outbursts of U Gem-type and SU UMa-type DN is A\ 
0.-18 + 0.03 with dispersion rms = 0.57. 

As discussed in iKotko et all (120121) . PTF1J0719 is the only 
system in the AM CVn class of binaries where short outbursts 
can be firmly classified as "normal" and the system considered as 
a helium counterpart of an SU UMa-type DN. With PTF1J0719 
taken into account, the coefficient of the linear fit remains almost 
unchanged. 

We conclude that the r de c - Pd.max relation is universal for 
normal outbursts of all classes of cataclysmic variables. To esti- 
mate the ah parameter one needs to find the relevant Td ec -Pd,max 
relation for model light-curves calculated with different and 
compare the result with observations. We chose four values of 
ah- 0.05, 0.1, 0.2, and 0.3, for each of them the set of models 
with different mass transfer rates M tr , primary masses Mi and 
maximum disc radii Rd.max were calculated. The decay rates of 
the synthetic outbursts were measured in the same manner as in 
the observational case. 
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To obt ain the synthetic light -curves we used the code de- 
scribed in IHameurv et al.l (Il998l) . which differs from the code 
used by S99; for example the input parameters are different^ 
and in our code the adaptive grid enables high resolution of the 
fronts. Such differences should not affect modelling of the decay 
phase of the outburst. 
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Rd,max[1.0x10 10 cm] 

Figure 1. The Tdec-^d,max relation. Filled symbols: U Gem-type systems 
(circles), normal outbursts of SU UMa-type systems (diamonds), and 
PTF1J0719 (square). Open symbols: models with a h = 0.1 (circles), 
ah = 0.2 (squares), = 0.3 (triangles). The dotted line - linear fit to 
the observational data (filled symbols) in the form T dcc = AiR imiix , with 
Ai = 0.48. (For the sake of clarity models with = 0.05 were not 
plotted.) 



When measuring Td ec of the outbursts for a wide range of 
model parameters, one has to pay attention to several problems: 

(a) For large discs (/?d,max > 5.0 x 10 10 cm) and high primary 
masses (M\ > 1 M ), so-called "refiares" appear during 
the decline from maximum. They are an indication of the 
cooling and heating front reflections in the disc, where the 
surf ace density is close to its critical valu e (for details see 
e.g. iMenou et al.ll2000h bubus et al.ll2001l) . In this case the 
model outburst cannot be considered as normal. 

(b) In large discs and for high values of ah , inside-out heating 
fronts may not be able to propagate up to the outer disc edge 
because such values of the viscosity parameter decrease the 
value of El. (see the formulae in AppendixlAl and, with in- 
creasing R, increase the possibility of a cooling front form- 
ing right behind the heating front. In this case the cooling 
front will start at R <K Rd.max and the decay rate will not be 
connected with the actual radial extent of the disc. 

(c) Models with the same parameters except for have differ- 
ent stability limits. 

The above-mentioned DIM properties have their reflection 
in the distribution of the model points in Fig. \T\ For large /?d,max 
models with higher tend to deviate more from the empirical 
Tdec-^d,max relation. 

2 S99 uses Mi and Mi (which define the orbital period of a semi- 
detached binary syste m) as the input parame ters and the mean disc ra- 
dius (Ra), while in the Ha meurv et al.l dl998h code, the input is Mi and 
(Rj). Smak defines only one value of a, while we take as the input pa- 
rameters both ah and a c . 



Linear fits to the Td ec = A a R^ md x relation obtained for the 
models with different aj, give 

1. Aoos = 1-624 ± 0.235 for or h = 0.05, 

2. A .i = 0.525 ± 0.128 fora h = 0.1, 

3. A 2 = 0.338 ± 0.036 for a h = 0.2, 

4. A .3 = 0.151 ± 0.031 foro; h = 0.3. 

The coefficients A a show a clear tendency to decrease when a 
higher is set in the model. The comparison with A\ obtained 
from the fits to empirical data shows that a^e [0.1, 0.2], with no 
unambiguous preference for one of these values, thus confirming 
conclusions obtained by S99. 

SU UMa-type stars and their superoutbursts provide another 
piece of interesting information. The decay from superoutburst 
may be divided into at least two phases - the plateau and fast 
decay phases. According to the enhanced mas s-transfer (EMT) 
model dKotko et aLlboTl ISmakll2008L I2009allbll3ldh . during the 
plateau phase the slow decline of the system luminosity is caused 
by accretion-driven depletion of the excess matter provided by 
the enhanced mass transfer from the secondary. This phase ends 
when a cooling front forms, and so the following fast decline is 
caused by the mechanism producing normal outbursts. 

Based on this, we measured Td ec during the fast decay phase 
of SU UMa superoutbursts and found that they are approxi- 
mately the same as the Td ec measured for their normal outbursts. 
The same is true of superoutbursts and normal outbursts in the 
models calculated with the pre scription for the M t , enhancement 
given in Ham eurv et al.l (1 19971) . The decay time in the fast de- 
cline phase was measured as the time interval between the time 
the system luminosity was 1 mag below the start of the decline 
phase to the time when the system was 2 mag below it. 

This conclusion is very promising in the context of evalu- 
ating ffh in AM CVn stars. As already mentioned, the normal 
outbursts in AM CVn stars are rarely detected, and the outburst 
cycle is dominated by superoutbursts. However, with well ob- 
served, fast decay phases of the superoutbursts in AM CVn stars, 
it will be possible to estimate in helium-dominated discs more 
precisely. Unfortunately the currently available data are not of 
sufficient quality to permit such investigations. 

3.2. The outburst width - orbital period relation 

van Paradii si (1 19831) showed that there exists a positive correla- 
tion between the outburst width W and the orbital period but 
concluded that narrow and wide outbursts should be considered 
separately. 

To consistently determine the outburst width in various sys- 
tems, the magnit ude level at which i t is measured has to be de- 
fined. Following Ivan Paradii si (1 19831) S99 defines W as the time 
interval during which the system luminosity is above the level 
set at 2 mag below th e outburst maximum . 

Using data from [van Paradiii dl983l) . S99 finds the coeffi- 
cients for 

- the linear dependence in the form W = Cw.swPmb'- Cty,s99 = 
1.39 ±0.06; 

- the non-linear dependence in the form W = Cw^Porb^ 99 '- 
C tv,S99 = 2.01 ± 0.29 and #399 = 0.78 ±0.11. 

Since, as in S99, we find the linear fit to be of superior quality, 
in what follows we do not use the non-linear fitting formula. We 
also only use narrow outbursts as "generic" normal outbursts. 
The linear fit to our data (18 systems) gives Cw = 0.99 ±0.12, 
so the agreement with S99 is not as good as for the decay times. 
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The analogous procedure applied to different subsets of our 
data gives 



- for U Gem stars from lAk et all (|2002|) only: C w = 0.79 ± 
0.07; 

- for U Gem stars from lAk et alj d2002l) complemented with 
our measurements: Cw = 0.90 ± 0.10. 

One concludes that Cw depends on the choice of the DN sample. 

The measurements of the outbursts width are clearly more 
vulnerable to uncertainties that are connected with the precise 
determination of the outburst maximum and with usually sparser 
data coverage of the outburst rise in comparison to the outburst 
decline. Moreover, except for systems observed intensively for 
a long time (such as the already mentioned SS Cyg or U Gem), 
straightforward assessment of which outbursts are narrow and 
which are wide may be problematic. Nevertheless, it is worth 
comparing observations with models as has been done in Section 
13. ll for Tdec ■ The width W of synthetic outbursts has been defined 
in the same manner as in the observational case. 

When applying the linear dependence W = C a R^ max to the 
same set of models as in Section[3~Tl one obtains 



1. Co.os = 3.222 ± 1.132 forcr h = 0.05, 

2. Co i = 1.794 ± 0.285 for a h = 0.1, 

3. Co 2 = 1.502 + 0.1 12 for a h = 0.2, 

4. Co.3 = 0.872 + 0.165 for a h = 0.3. 

Comparison of the above listed C a 's with C\ = 1.496 + 0.135 de- 
termined from the observed relation W = C\R^ max again favours 
a^e [0.1, 0.2] with even stronger indication of a^- 0.2. 

The model and observational data with fitted linear depen- 
dence between W and /?d,max are presented in Fig. [2] Even if 




4 5 6 7 

FW<[1-Ox10 10 cm] 

Figure 2. The W-Ri mm relation. Filled symbols: U Gem-type systems 
(circles), normal outbursts of SU UMa-type systems (diamonds), and 
PTF1J0719 (square). Open symbols: models with a h = 0.1 (circles), 
ah = 0.2 (squares), aj, = 0.3 (triangles). The dotted line - linear fit to 
the observational data (filled symbols) with a coefficient Ci = 1.496. 
(For the sake of clarity models with = 0.05 were not plotted.) 



W is not determined well enough to provide a firm value of ah 
from the W-Rd,max relation, the results obtained totally preclude 
a h <sc 0.1. 



3.3. The decay time from the DIM 

Since we have assumed that the decay from outburst's maxi- 
mum is described by the standard DIM, it is worth checking how 
strongly model-dependent the results are concerning the value of 
the viscosity parameter. In principle, the setting is very simple: 
shrinking (from outside) hot-disc configurations decay through 
a set a quasi-stationary solutions corresponding to monotoni- 
cally diminishing accretion rates. This prompted S99 to estimate 
the decay time as: 



R, 



d,max 



(3) 



where the viscous speed v v i s ~ v/R. However, S99 used 
a non-standard definition of the kinematic viscosity coeffi- 
cient. To clarify thi s poin t we recall that the idea behind the 
Shakura & Sunvaevl d 1973b a ansatz is that the component r r< p 
of the stress-tensor should be proportional to the total pressure 
P with the proportionality constant a describing the efficiency 
of the angular momentum transport due to turbulence in the 
disc: t>0 = -aP, where < a < 1. On the other hand, in 
a differentially rotating fluid the tangential stress is defined as 
7>0 = r/dQ/dR, where 77 is the dynamical viscosity coefficient, R 
the radius, and Q. the angular speed. For a Keplerian disc, there- 
fore, the kinematic viscosity coefficient v = rj/p, where p is the 
density, can be written as 



2 c z s 
v - —a — — 

3 CIk 



2 

—aCoH 

3 ' 



(4) 



where c s is the sound speed, the Keplerian angular speed., 
and H * c s /Q.k the disc (semi) height-scale and not the disc's 
actual (e.g. photospheric) height zo as assumed in S99. 

Fortunately, was not calculated in S99 directly from the 
analytical formula but obtained (as we did in the previous sec- 
tions) from comparing the theoretical and empirical dependence 
between the outburst decay rate (or the outburst width) and the 
orbital period, so this article's conclusions concerning the value 
of the viscosity parameter in hot accretion discs remain valid. 
Nevertheless, it would be useful to derive an analytical formula 
for rdec, which would only depend on and observables since it 
would allow dispensing with using numerically calculated mod- 
els, especially because these must be selected according to cri- 
teria discussed in Sect. 13. II In Table|2]we compared decay rates 
calculated from numerical models (Td ec m ) with the values (Td ec v ) 
corresponding to Eq. 0). (In the selected models fronts always 
propagate through the whole disc extent.) Clearly the rates ob- 
tained from Eq. (O are much too long and the approximation 
used was too crude. 

Indeed, the decrease in luminosity after the outburst maxi- 
mum is the effect of two mechanisms: (1) the depletion of the 
matter from the disc due to (viscous) accretion onto the central 
object, and (2) the propagation of the cooli ng front through the 
disc (also a viscous process, for details see Men ou et al. 1999; 
Lasota 1200 lh . Taking the disc shrinking into account lowers the 
decay time. 

W assume that the decay time fd ec is the time it takes the 
system luminosity to drop from the maximum to the quiescence 
level and that it may be written approximately as 



fde 



Ra 



|Vde. 



where (as before) R^ yI 
decline velocity. 



(5) 



is the maximum disc radius and Vd ec a 
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Table 2. Decay rates for several solar composition disc models. 



Model 


a h 


Mi 
(M ) 


"d.max 
/into 

(1U cm) 


An 

(mag) 


Tdec™ 

(d/mag) 


_ i 

Tdec 

(d/mag) 


Tdec'' 

(d/mag) 


1. 


0.2 


1.0 


4.18 


6.3 


1.62 


1.17 


8.38 


2. 


0.2 


1.0 


1.29 


4.1 


0.78 


0.6 


7.17 


3. 


0.2 


0.6 


1.29 


5.6 


0.48 


0.42 


7.65 


4. 


0.2 


0.6 


0.76 


6.3 


0.23 


0.23 


2.75 


5. 


0.3 


1.3 


3.04 


3.3 


0.97 


1.13 


10.37 


6. 


0.3 


1.0 


1.29 


3.2 


0.52 


0.51 


6.1 


7. 


0.3 


1.0 


0.78 


3.7 


0.29 


0.27 


4.4 


8. 


0.3 


0.6 


0.77 


3.7 


0.16 


0.26 


3.14 



Notes. Tdec'" is measured from the model, Tdec' calculated from Eq. ifTJ, and Tde C ' from Eq. <[3j . 



To find a simple formula we assumed that Vd ec is the super- 
position of (a) the inward viscous velocity v v ; sc of matter in the 
hot part of the disc into which the cooling front pr opagates and 
(b) th e outward velocity of the gas at the front V/J_| (Menou et al. 
1999) allowing the inward propagation of the front; i.e. |vd ec | ~ 
|v v isc + Vf|. This crude approximation treating the two veloci- 
ties as an average over space values gives quite good results, as 
seen in Table [2] This is b ecause for most of the disc's extent (in 
the "asymptotic regime": iMenou et al.ll 19991 their Fig. 7), Vp is 
roughly constant. Since v v i sc « v/R, from Eq. (0 one has 



(6) 



where P is the total pressure and c s - ^/P/p. 

From numerical simulations, Menou et al. (1999) find that 
Vf ~ 1 /7ar n c s . We confirmed this result for cooling front veloc- 
ity in solar composition discs and found that it also applies to a 
helium-dominated disc. The final formula for fd ec is thus: 



1R A 



a hCj (l + 14/3(c s /v^(« d , max )))' 



(7) 



The speed of sound can be expressed in terms of the central tem- 
perature in the disc: c s = y/kT c /mH, where k is the Boltzmann 
constant and mn the hydrogen molecular mass. (In the case of 
a helium disc it should be replaced by helium molecular mass). 
The numerical fit to the temperature at the cooling front found 
from models of solar-composition discs gives 



T c *4.7 x 10 4 K, 



(8) 



with no dependence on disc parameters. From Table [2] one can 
see that Eq. can give reliable estimates of of the viscosity 
parameter, certainly better than Eq. (0. 

For AM CVn stars, e.g. PTF1J0719, the chemical composi- 
tion (Y = 0.98 Z = 0.02) gives 



1.1 x 10 3 K 



(9) 



From Eqs. and (0 (or (0) it is clear that ?d ec depends on 
^Amax, Mi , and ah . The primary mass Mi de termines the whit e 
dwarf radius R\ through the M - R relation dNauenberglU"972l) . 



3 We follow here the Menou et al.l (1999) notation according to which 
v f is not the front velocity but the gas velocity at the front. The ratio 
of the gas velocity at the cooling front to the cooling front velocity is 
typically ~ 2 



Both Mi and ^d.max define the disc's extent since the model as- 
sumes that the inner disc radius R m — R\. 

To compare the observed outbursts decay rates Td ec ° with the 
analytical decay rates Td ec ', the derived fdec has to be divided by 
the amplitude of the outburst. For each system with measured 
f or b and estimated Mi and M2 (necessary for calculating ^d.maxX 
there is only one free parameter left: a^. Thus the conformity 
between the observed decay rate Td ec ° and Td ec ' calculated from 
Eq. (0 can be attained by adjusting ah. We assumed = 0.2. 
The calculation results are summarized in TableQ] It is seen that 
the calculated time is close to the observed one in most cases. 
In a few cases the discrepancy is large. It is not clear if it is due 
to the imprecision in measuring the decline time, the peculiar 
nature of the outbursts, or to the non-universal value of a^. After 
all, since MRI does not give the correct value of this parameter, 
we do not really know what physical mechanism drives accretion 
in hot dwarf nova discs, so there is no reason to assume that it is 
"generic". 

One should keep in mind, however, that the significant im- 
pact on a n has the determination accuracy of the observed out- 
bursts amplitudes. An underestimate of A n may be the cause of 
significantly higher a h for TW Vir, WX Hyi, and PTF1J0719. 



4. The Kukarkin-Parenago relation 

The first to suggest a relation between the outburst amplitude A n 
and th e outburst recurrence time T n were iKukarkin & Pare nago 
(1934). However, since their sample contained outbursts of 
both re current and dwarf novae, its reality has been ques - 
tioned (iPavne-G aposchkrnl 1 1 95 7l 119771 iBath & Shavivl ll978). 
In his seminal revi ew article, Smakl ( 1984b ) attributes the cor- 
relation to |Pavne-"G aposchkin) (119771) and stresses its statistical 
character. Finally, a nalysing dwarf nova normal-outbursts data, 
Ivan Paradiisl (IT985) concluded that 'the amplitudes and average 
recurrence times of dwarf novae are correlated". The most re- 
cent versio n of the Kukarkin-Parenago relation (hereafter K-P 
relation) in Warner (2003) takes the form 



A n = (0.7 ±0.43) + (1.9 + 0.22) log T a , 



(10) 



where T n is in days and A n in magnitudes. 

Although it has been argued that the K-P relation might 
represent some glo bal and average properties of DN outbursts 
(Ivan Paradii sll 1985b . to the best of our knowledge no derivation 
from the model has been attempted until now. 
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While the other relations, connecting various quantities 
characterizing the binary systems and their outburst light- 
curves (such as the absolute visual mag nitude-at-max imum 
Afym axCP orb) or the TdecC^orb) relations (see I Warner! I2003L and 
Sect. 13.31 ) follow directly from the DIM, derivation of the K-P 
relation is not straightforward. Of course as the outburst ampli- 
tude is related to the mass of the disc and the mass of the disc 
to the accumulation time, a relation of the K-P form sh ould be 
expected in principle. For example. Ivan Paradjjl d 1985b specu- 
lated that the average amount of mass AM qu j es transferred during 



an average recurrence time T n : AM, 



quies 



r n M, 



quies 



is constant 



over the dwarf nova population, but also warned about possible 
selection effects and model dependence. Here, we try to examine 
what kind A n (T n ) relation, if any, can be deduced from the DIM 
in its simplest (and when possible simplified) form. 

For simplicity we assume that, during quiescence, the accu- 
mulation rate M accum of the mass in the disc is approximately 
equal to the mass transfer rate from the secondary M tl (this as- 
sumes no truncation of the inner edge of the disc; no "leaky" 
disc) and that the mean accretion rate during outburst is about 
half the maximum accretion rate during outburst (M out b) ~ 
5^accr,max- (This follows from the shape of the function M accr (f) 
in the model, which is rapidly rising and then approximately ex- 
ponentially decreasing during the outburst.) 

It is worth remarking here that the maximum accretion rate 
A^accr.max is only approximately equal to M+ it (^d,max) ar, d sim- 
ulations clearly show that the mass accretion reaches the max- 
imum only after the cooling front has started to propagate (see 
Fig. [3]). The accretion rate keeps rising when the heating front 
arrives at the outer disc edge and a cooling front starts propagat- 
ing. This is because just before the launch of the cooling front, 
the E-profile had not yet reached the stationary hot disc £ shape 
and the mass from the heated parts near the outer edge had not 
have time to diffuse fully inwards. During the initial phase of 
the cooling-front propagation, this mass excess keeps diffusing 
inwards since it is also "shoved" by the incoming cooling front. 
The mass accretion rate will eventually drop due to the appear- 
ance of the mass shortage in the inner parts of the disc caused 
by two mechanisms: (1) accretion onto the white dwarf and (2) 
the strong outflow of mass at the cooling front, which shuffles 
the mass to the outer parts of the disc (see Sect. 13.31 ). Guided 
by the simulation results, we assume M accr max » eM^ (/?d,max)> 
with e ~ 3 in the following. 

The amount of mass accreted during the outburst decay is 
equal to the mass accumulated in the disc during quiescence and 
the rise to outburst maximum: AM accr = AM accum . 



(^outb)fdec — ^tr(fquiesc + frfse) 



(ID 



where f qu j e sc is the duration of the quiescence, f^ the time it 
takes the outburst to reach its (bolometric) luminosity maximum, 
and fdec the duration of the outburst decay to the quiescence level. 
Strictly speaking, mass accumulation also occurs during the last 
part of the decline from maximum (from points 1 to 2 in Fig. |4j 
but this has a negligible effect on the total mass balance. 

The decay time fdec is calculated as described in Section [3~3l 
The outburst recurrence time T n is counted from the onset 
of the outburst to the onset of the following one (T n = fn Se + 
fdec + fquiescX while the accumulation time is f qu iesc + frfse- fee 
is non negligible in some type of outbursts, and it is important 
to include it to be able to account for e.g. the "cycling state" 
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Figure 3. Top: Evolution of the £ profile during the propagation of 
the cooling front. Numbers next to the lines stand for the subsequent 
moments of the front propagation (1 is near the moment of the inset of 
the cooling front). Bottom: Evolution of the M accr profile. The numbers 
correspond to the same time points as those in S. 



outbursts where no quiescence phase is present.) Substituting T a 
and f accum to Eq. CO} gives 



T„ = 



(Mouft) 

M tr 



+ 1 fdec = 



M a< 



2M a 



+ 1 fd, 



(12) 



Since M accrmax * e -^Cit ^ d -* an< ^ me instability condi- 
tion requires Mtr < MX (Ri), the ratio in the brackets is 
M+ it (R d ) /(2M tr ) » 1 and one can take 



2M tr 



(13) 



(in what follows we drop the index "max" in Rj). 

Assuming that Ri/Rd <Z 1 , the luminosity a t outburst maxi- 
mum can be approximated as (Frank et al. 2002) 



GM x eM + mi {R A ) 



(14) 



where R\ is a white dwarf radius and Mi its mass. The luminos- 
ity at minimum light can be estimated from the model as 



2R A ~g 



(15) 
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Figure 4. The changes of the disc mass M iisc (top), magnitude (mid- 
dle), and mass accretion rate M accr at Ri n (bottom) during one out- 
burst cycle. The horizontal red line at the M xa (t) plot is the constant 
mass transfer rate M tr , and the horizontal green line is M^. it (/? dmax ) for 
this model. The parameters of the model are = 0.2, a c = 0.05, 
Mi = 1.0 Mo, M u = 6.0 x 10 16 g/s, and (R d ) = 1.2 x 10 10 cm. The 
point number 1 on all three plots refers to the point where the mass ac- 
cumulation in the disc starts, point 2 is where the mass depletion from 
the disc starts, and point 3 is the end of this cycle when the mass starts 
to accumulate in the disc again. 

Table 3. Six sets of parameters for which the theoretical K-P 
relation has been calculated. 



Model 


Ci 


ah 


"d,max 


Mi 


1 


-2.5 


0.01 


2 


1 


2 


-0.95 


0.2 


0.6 


0.6 


3 


-0.2 


0.2 


2 


1 


4* 


0.1 


0.2 


0.7 


1 


5 


1.2 


0.2 


5 


1.2 


6* 


1.3 


0.2 


0.8 


1.2 



Notes, ov hot disc viscosity parameter, and C\ - the constant from Eq. 
d!7t calculated for a given set of parameters. Models 2 and 5 correspond 
to the lower and upper limits of the theoretical K-P relation for solar 
discs with a h = 0.2; "*" models are calculated for helium discs (Y = 
0.98 Z = 0.02): model 4* gives the K-P relation for PTF1J0719, model 
6* gives the K-P relation for CR Boo and V803 Cen in their cycling 
states. The models are plotted with lines in Fig. [5] 



where g~2 (lldan et al.lll999h . 

The amplitude A n is the difference between the magnitudes 
at maximum Mv. mm and minimum My.min given by 

A„ = M v ,_ - M v , max + BC1 = 2.5 log + BC + _ (16) 

where BCt is the difference between the bolometric corrections 
at maximum and minimum. From Eqs. ( [TBI . ( fT4b and ( fl3b we 
thus get 

A„*C, +2.5 log T a (17) 
where 

Ci = 2.5 log 2g - 2.5 log f dec + BC max - BC mm . (18) 
The last step is to estimate the bolometric corrections BC max and 

BC Ynin- 

Fro m the definition, BC max = M DO i,max - My, max . Following 

Smak| dl989l) we calculate the visual luminosity at maximum 
£v,max with spectral energy distribution of a black body disc in- 
tegrated over the visual band frequencies. The luminosity Lv.max 
(and so BC max ) depends on Mi, R\, Ra and M accr , max . The bolo- 
metric correction BC m i„ can be estimated from the spectral en- 
ergy distribution calculated for quiescent disc models with ef- 
fective temperature < 5000 K. We us ed spectra calc ulated by 
Irit Idan (private communication; see Idan et al.ll2010h . Based 
on that we assumed BC 1Tun « -0.4. 

Equation $T7\ does not correspond exactly to the K-P rela- 
tion (see Eqs. [T0l[T9l l. The slope of the theoretical relation is al- 
ways 2.5 (by construction) as compared with ~ 2 obtained from 
fits to observations. Considering the typical spread of parame- 
ters, C\ is contained between ~ -1.0 and ~ 1.2 in hydrogen- 
dominated discs, and is between 0.2 and 1.3 in helium discs 

(see Tabled, to be compared with C\ between 0.3 and 2.0 for the 
K-P relation. Considering the very large scatter of observational 
data, this can be considered a fairly satisfactory result. This is 
only true of the theoretical relation obtained assuming = 0.2. 
The same relation with = 0.01 gives a totally unacceptable 
representation of the A n (T n ) relation, thus confirming the con- 
clusion of the previous sections that aj, ~ 0.2. (Compare models 
1 and 3 to see how the K-P relation changes when is decreased 
from ofh = 0.2 to = 0.01 in a solar disc.) The various A n (T n ) re- 
lations are plotted in Fig. [5] We also marked on this diagram 
the values of A a (T n ) for a subset of dwarf novae and outburst- 
ing AM Cn stars. The sample of the systems p resented in Fig . 
[5] consists of U Gem-type binaries taken from lAk et alj (2 002) 
and li sted in Table[T] SU UMa-type binaries taken from lAk et aTl 
(2002) (also liste d in Table[D and from the updated Cataclysmic 
Binaries Catalog dRitter & Kolbll2003h and three AM CVn-type 
stars for which A n and T n were measured from their light-curves. 
For CR Boo and V803 Cen th e measured A n and T n relate to the 
outbursts in the cycling state (Patt erson et al.ll2000l) . 

In our sam ple 15 systems out of 43 are the same as used by 
IWarnerl (2003). The linear fit to our sample gives 

A n = (1.3 ±0.6) + (1.6 ±0.3) log T,, (19) 

Upper and lower uncertainties of this relation are marked in 
FigS 

Some of the U Gem-type binaries in Fig. [5] are marked with 
asterisks and their names as examples to show how observed 
systems correspond to the theoretical lines. For the same purpose 
the name of one of AM CVn stars is shown on the plot. 

The independence of the theoretical K-P relation from M tr 
is the consequence of the assumption M accr max ~ eMij t (7?d,max) 
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Figure 5. Various A n (r n ) relations. (Numbers in brackets above the 
lines correspond to the row numbers in Table [3j- The solid lines (2) 
and (5) correspond respectively to the "lower and upper limits" (see 
text) deduced from the theoretical relation for solar-aboundance discs 
with a h = 0.2 (Ci = -0.95 and 1.2). Helium disc theoretical rela- 
tions are represented by the dotted lines (4) and (6) (Ci =0.1 and 1.3). 
Solid lines 1 and 3 illustrate the a h dependence of the A n (T n ) rela- 
tion: for the same binary parameters, they correspond to = 0.01 
and ff h = 0.2, respectively. The upper and lower uncertainty of ob- 
servational K-P relation fitted to the systems are marked with thick, 
dashed lines. The sample of binaries marked on the plot consists of 
U Gem-type systems listed in Table [T] (full circles and asteriks), AM 
CVn systems: PTF1J0719, CR Boo and V803 Cen (both in cycling 
state outbursts) {op en circles ) and SU UMa stars (normal outbursts only, 
/;/// squa res) fro m Ak et al. (2002) and updated Cataclysmic Binaries 
Catalog dRitter & Kofbll2003l) . 

since Mt^XR d imax ) does not depend on M a . Systems with larger 
(more extended) discs and more massive Mi have higher A„ for a 
given T n than systems with small discs or less massive primaries 
(compare lines (2) - (5), (4) - (6)). 

Despite the simplifications and approximations assumed in 
the derivation, the theoretical K-P relation follows the observa- 
tional data reasonably well. One concludes that normal dwarf 
nova outbursts are indeed the results of filling and emptying of 
an accretion disc, as assumed in the model. The parameter that 
has the deciding influence on the recurrence time and amplitude 
of normal outbursts is the disc's extent. 

5. Conclusions 

As in S99, the main conclusion of the present paper is that 
in ionized dwarf-nova accretion discs the viscosity parameter 
ffh ~ 0.2. The same conclusion is presumably also true for 
helium-dominated discs in outbursting AM CVn stars, although 
there the statistics on which it is based are fairly poor. Although 
there is no evidence that the value of «h is universal, it can be 
firmly established that, even if it varies over the cataclysmic vari- 
able population, it cannot be as low as 0.01, the value resulting 
from numerical simulations of the MRJ, which is the mecha- 
nism that is supposed to drive accretion in hot Keplerian discs. 
Therefore it is not preposterous to suggest that solving this dis- 
crepancy between observations and theory should become the 
main subject of interest of researchers studying disc accretion 
mechanism^. 

4 After the submission of th e present article, a paper on this subject 
bv lLatter"&^ Papaloizou (2012) has been posted on astro-ph. 
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Appendix A: Critical parameters for solar 
composition discs. 

The method of obtaining the formulae and the formulae themselves can be found 
in lLasota et ail <2008l) . 

For the solar composition hydrogen, helium and metal mass fractions X = 
0.7 Y = 0.28 Z = 0.02 one obtains from fits to S-curves 

S + = 39.9 ao ( j- 80 R; o 11 m7 - 37 gcm- 2 



9 



I. Kotko & J.-P. Lasota: The viscosity parameter a 



Z- = 74.6 o^" 83 /4 18 m7°- 40 gcrrr 2 

T+ = 30000 a ( J 18 K? 04 m7 01 K 

T- = 8249 a° \ 4 R^ A m? 04 K (A.l) 

T+ ff = 6890/? w 00,) m»- 03 K 

T eff = 521()« ro 01 mrK 

J>+ , m „inl5 ,-0.01 D 2.64 -0.89„ o -l 
M = 8.07x10 a Ql K| njj gs 

Jl>- T A/1 „0.01 r,2.58 „,-0.85„ -1 

M = 2.64 X 10 ci^ 1 m 1 gs , 

where are the critical surface densities for the hot (+) and cold (-) thermal- 
equilibrium solutions. Similarly, r*, T* ff are the critical values of mid-plane and 
effective temperatures, while M± correspond to critical accretion rates, m\ is the 
primary mass in solar units and i?io is the radius in units of 10 10 cm. 
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